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Part I 


INTRODUCTION 

The mass spectrometer has long been used as an analyt- 
ical tool in many disciplines. Relative to the space program, 
the mass spectrometer was one of the first major analytical 
instruments to be utilized from the very inception of the 
nation's program. 

Virtually every major rocket was equipped with a mass 
spectrometer to measure the atmosphere either on ascent or 
descent and in all recent missions the mass spectrometer has 
provided a wealth of information which cannot be obtained with 
any other type of instrumentation. Many types of mass spec- 
trometers have monitored the upper atmosphere of the earth 
and are continuing to do so. Additionally, deep space probes 
can be expected to continue to utilize mass spectrometry, and 
the surfaces of planets are amenable to analysis by mass 
spectrometry via telemetering of their chemical composition. 

Another major area in which mass spectrometry has made a 
contribution to NASA sponsored programs relates to the high 
purity materials that have been developed for solid state 
electronics. Without mass spectrometr ic analysis to develop 
a solid state technology, the rapid advance of transistors 
and other microcircuitry, computers, etc. would not have 
been possible, inasmuch as the mass spectrometer has a sensi- 
tivity for analyzing trace elements that generally exceeds all 



other analytical methodologies. 

The present work, however, relates to a more generic 
type of materials research that has applicability to high 
strength metals and alloys for advanced space vehicles, and 
even to the ultimate performance of composites in supersonic 
transports. Hence, this research represents an even further 
extension of the manifold applications for mass spectrometry 
that have occurred during the last two decades. 

Specifically, this paper addresses the problem of the 
high temperature diffusion in metals, and it does so in a 
very unique manner. Until the last few years, there was 
rather little interest in the trace elements of metals as 
opposed to some of the major constituents that are generally 
recognized as alloying elements. However, we are now approach- 
ing an era in which the efficiency and reliability of many 
systems, e.g., freedom from corrosion, the attainment 
of theoretical yield strengths, and the macroscopic properties 
of materials generally are dependent upon trace metals and 
their transport to grain boundaries. Hence, while this present 
research effort is focused on a few specific materials, the 
technology that it utilizes should have a very broad scale 
application, and it is in this context that we are presenting 
results of a preliminary nature which will hopefully lead to 
an extended future study of diffusion in many engineering 
materials . 

We are pleased to acknowledge the technical input and 
general support from Dr. George M. Wood of the NASA-Langley 
Research Center during the entire period of this investigation. 
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and for research funding from the National Aeronautics and 
Space Administration under Research Grant NSG-1360. 

The first basic theory for calculating diffusion coeffi- 
cients was presented by Adolf Pick in 1855.^ Pick who modified 

2 

Pourier's heat conduction equations, hypothesized that in 
an isotropic medium the quantity of diffusing substance Q, 
which passes in a unit of time through a unit of transverse 
cross sectional area, is proportional to the concentration 
gradient measured along the normal to the section: 

Q = -DCdC/dx) . Cl) 

The above equation is called pick's first law for equilibrium 

flow, C is the concentration of the diffusion substance in 

units of mass/cm^, x is the spatial coordinate in units of 

2 

cm., and D is the diffusion coefficient with units of cm /sec. 
In the case of non- equilibrium flow. Pick's second law can be 
derived from the first law by considering the rate of accumula- 
tion of the diffusing substance in a given element of volume 
as the difference between the incoming and outgoing fluxes 
per unit time. This yields 


^ = D ^ * D ^ - D ^ 

dx^ dy'^ dz^ 


C2) 


Pick's laws were shown to be approximations to more 
general transport equations, and they may apply only to several 
atomic distances from solute sources or sinks, when the solvent 
is a homogeneous medium and when the solute concentrations are 

3 

small. Boltzmann extended the application of Pick's equations 
to large solute concentrations,^ and these modified equations 
were graphically solved by Matano^ in 1936 for metallic 
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diffusion. To date, experimental determinations of the diffu- 
sion coefficient have always followed a prescribed procedure. 
The solute material was allowed to diffuse at a known temper- 
ature and for a known time into a block of solvent material 
of such a shape that Pick’s laws yielded an exact mathematical 
solution for the diffusion process in the system. Next, the 
solute concentration was measured at known points throughout 
the solvent. These values of solute concentration, and the 
diffusion coefficient was found algebraically or graphically. 
The problem of experimentally determining the diffusion co- 
efficient thus became one of accurately determining the solute 
concentrations throughout the solvent, and of fulfilling the 
initial and boundary conditions applicable to a particular form 
of Pick’s laws. 

The methods which various experimenters have used to 
determine the solute concentrations and to satisfy the Pick’s 
law boundary conditions in their experiments are quite diverse. 
An excellent summary of these different methods have been 
compiled by Gertsriken and Dekhtyar.^ 

The introduction of radioisotopes in the 194Q’s enabled 
the concentration of a radioactive solute in a solvent 
specimen to be determined quite readily by nuclear radiation 
detectors. It is now advantageous to illustrate these 
standard methods. 

Serial Sectioning ; This method involving radioactive tracers 
is most frequently used in the study of diffusion. The sample 
after being held at a given temperature for a known amount of 
time (annealing time) to allow the radioactive tracer to 
diffuse into the sample, is subjected to a series o*f sectioning 
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operations. The sample is weighed before and after each 
cut to determine the thickness of the material removed. 

The cuttings are collected and their activity measured with 
a Geiger or scintillation counter. From the weights and 
the density, the coordinates of the midpoints of each slice, 

(x) , are calculated. A plot of the natural logarithm of 
the activity versus x has a slope of -l/4Dt for volume 
diffusion. Hence, if the time of the anneal, t, is known, 
the diffusion coefficient D may be calculated. 

The main advantage of the serial sectioning method is 
that it is simple, direct, and doesn’t depend on the properties 
of the radiation from the radioactive material used, assuming 
no damage to substrate by ionizing radiation from the radio- 
isotope . 

Residual Activity Method : This method is similar to the 

serial sectioning analysis in that layers of thickness x 
are removed, but the total remaining activity, I, of the 
sample is measured in this case. 

In the case of weakly absorbed radiation (strong gamma 

2 

rays) a plot of ln(-dl/dx) versus x should be a straight 
line of slope -l/4Dt; and in the case of strongly absorbed 
radiation (weak beta-rays} a plot of ln(I) versus x^ should 
also be a straight line of slope -l/4Dt. 

Surface Decrease Method : In this method, the total activity 

of the specimen is measured as a function of time. No section- 
ing of the sample is necessary. 
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Autoradiography : The principle of this method is to determine 

the photographic density of the blackening of an exposed 
emulsion as a function of the distance from the interface of 
the specimen with the initial radioactive deposit. The speci- 
men is usually cut at a measured angle a (usually 90®) to the 
initial face, and the cut face is placed in contact with a 
piece of appropriate X-ray film. The blackening of the emulsion 
(which is usually found to be directly proportional to the 

cconcentration of activity) is measured with a micro-densitom- 

2 

eter. A plot of In (photographic density) versus x has a slope 
of -l/4Dt from which D may be calculated if the anneal time, 
t, is known. This method offered two distinct advantages; 1) 
since the filament was small, the times involved in heating 
and cooling it were small, and 2) because of the source mount- 
ing and the high vacua of the mass spectrometer the diffusing 
system was negligibly affected by any reaction between it and 
the surrounding environment at the experimental operating 
temperatures . 

With the advent of sophisticated electronics, new devices 
have been developed for studying solid state phenomena. The 
electron-beam microprobe is one of the most important as far 
as diffusion studies are concerned. The device collimates a 
beam of electrons into a 1 u diameter "pencil" which is 
directed onto the specimen surface at the spot to be analyzed. 
The electron bombardment causes characteristic X-rays to be 
mitted from the sample. These X-rays are then detected and 
analyzed to determine composition of the bombarded surface. 
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grains for elements 


The limit of detectability is about 10 

7 

of high atomic number. 

Existing procedures and techniques used for material 
characterization are adequate for materials which are to be 
operated in a normal temperature environment. However, at 
elevated temperatures (above 800°C) , there is a need for 
improved measuring techniques. It is very important to under- 
stand the effects of high temperature operation on the materials, 
i.e., to determine whether or not the elevated temperature 
will cause matter from the environment to be absorbed into the 
material and then be transported through it, or whether some 
constituent of the material will be removed. The transport of 
constituents either into or out of the material can change the 
structural integrity of the material. As the high temperatures 
diffusing system tends to react with its environment. Also 
rapidly heating and cooling of the solute-solvent system will 
often adversely affect the diffusion measurements. In 1960, 

g 

McCracken and Love developed a system for determining diffu- 
sion coefficients utilizing a mass spect rometer , and a thermal 
ionization source. 

Using a mass spectrometer with a thermal ionization source, 
only solute atoms of a single isotope would be observed, and 
diffusion rates in a specified temperature range could be ob- 
tained. The assumption is made that pick's solution for one 
dimensional diffusion in an outgasing thin plate can be applied. 
The solute ion current out of the filament as a function of 
time is thus : 
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■^ion ” ^ 


(3) 


where A is a constant of proportionality', D is the diffusion 
coefficient, x is the filament thickness, and t is the 
measurement time. Hence, the slope of a plot of the natural 
logarithm of the solute ion current versus time will yield 
the desired diffusion coefficient, provided the thickness of 

9 

the filament is known. 

General Advantages of the Mass Spectrometric Technique : There 

were many reasons for choosing the mass spectrometric approach 
for determining diffusion characteristics. As the mass 
spectrometer analyzes one mass at a time, there is no ambiguity 
caused by interfering masses. Since the filament is small, 
the times involved in heating and cooling it are also small. 

The probability of a reaction occuring between the high tem- 
perature specimen and the surrounding environment is negligible 
due to the mounting and the high vacua of the mass spectrometer. 
Finally, the sensitivity of the instrument permits the impurity 
level to be less than a part per billion, thus no doping of the 
sample is necessary and commercial grade material can be 
utilized . 

Importance of Sample Materials : The materials studied in this 

work were tantalum, zircaloy-2, and 304 stainless steel. These 
materials were chosen for two basic reasons. First, they 
represent a family of materials which go from being ultra pure 
(99.996% Ta) to a 1.5% alloy with zirconium, to a 70% Fe, 
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18.5% Cr, 9% Ni, alloy (304 stainless steel). Hopefully, 
this will enable one to draw conclusions about diffusion 
rates in alloys rather than ultra-pure single crystals. 
Secondly, these materials are of current interest to the 
nuclear industry where higher operating temperatures mean 
higher efficiencies. 

The alkali metals were chosen as the elements of inter- 
est primarily due to their relatively low ionization potential 
(making them more suitable for thermal ionization) and due 
to the fact that they probably appear in most materials as 
naturally occurring impurities - thus eliminating doping of 
the host matrix. It was decided that concentrating these 
experiments on an entire chemical family might yield an over- 
all diffusion mechanism which would allow diffusion rates of 
other impurities and sample materials to be approximated. 

Use of trade names or names of manufacturers in this 
report does not constitute official endorsement of such 
products or manufacturers, either expressed or implied, 
by the National Aeronautics and Space Administration. 
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Part II 


ATOMIC DIFFUSION THEORY 

A. General Considerations 

Diffusion is one of the mechanism by which matter 
is transported through matter. From the theory of specific 
heats, it is known that atoms in a crystal oscillate around 
their equilibrium position. Occasionally these oscillations 
become large enough to allow an atom to change sites. The 
net result of many such random movements of a large number of 
atoms is actual displacement of matter, the movement being 
activated by the thermal energy of the crystal. The path 
of an individual particle is an unpredictable zigzag. The 
length of an individual step in this zigzag is determined 
in a solid by its jump length, b, i.e., the distance the 
particle is able to move when it acquires enough energy to 
make one jump. 

The frequency of jumping, f, is the vibrational 
frequency multiplied by the probability that the vibrating 
atom has the activation energy to make the jump; 

f = (l/3)u) exp (Gq/RT) (4) 

where oj is the vibrational frequency, 1/3 allows roughly for 
the fact that only about 1/3 of the randomly directed vibra- 
tions are in the required direction, and G^ is the free 

11 

energy of activation. From this we can obtain an equation 
for the diffusion coefficient, D, 

D = (1/2) b^f (5) 
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B. Diffusion Mechanisms 


How the substance diffuses through the material is 
just as important as the diffusion coefficient formula itself. 
There are quite a variety of mechanisms by which an atom 
can move from one position to another in a crystalline 
structure . 

Vacancy Mechanism : in all other than perfect crys- 

tals, some of the lattice sites are unoccupied. These unoccu- 
pied sites are called vacancies. If one of the atoms on an 
adjacent site jumps into the vacancy, the atom is said to have 
diffused by a vacancy mechanism, see Figure 1(b). 

Interstitial Mechanism : When small atoms dissolve 

in metallic lattices as impurities so as to occupy inter- 
stitial positions between solvent atoms and an atom passes 
from one interstitial site to one of its nearest neighbor 
interstitial sites without permanently displacing any of 
the matrix atoms an interstitial mechanism is said to have 
taken place. The interstitial mechanism is thought to 
operate in alloys for those solute atoms which normally 
occupy interstitial positions. It will be dominant in any 

nonmetallic solid in which the diffusing interstitial 

1 2 

doesn't distort the lattice too much, see Figure 1 (a). 

Interchange or Exchange Mechanism : This is the 

direct exchange of two nearest neighbor atoms. 

Ring Mechanism : A more general form of the ex- 

change mechanism consisting of a number (three or more) of 
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atoms forming a closed ring. Atomic diffusion then takes 
place by rotation of the ring, 

Interst itialcy Mechanism : An interstitial atom 

moves from an interstitial site to an adjacent normal 
lattice site displacing the atom that had occupied that site 
into a new interstitial site, see Figure 1 (c) . 

Crowdion Mechanism : A line imperfection consist- 

ing of n nearest neighbor atoms compressed into a space 
normally occupied by atoms. Diffusion then takes 

place by movement along the line of atoms, see Figure 1 (d) . 

Dislocation Mechanism : Dislocations can provide 

easier paths for diffusion than a perfect lattice. The 
diffusion of dislocated atoms can be produced in two ways: 

1) the diffusion of the dislocated atoms over the inter- 
stices and 2) displacement of the dislocated atom from the 
interstice into a normal site of the lattice, and the atom 
situated in the lattice goes over into the interstice. 

C . Theoretical Methods of Diffusion Coefficient Calculations 

While the mechanisms of diffusion may be summarized 
from purely qualitative consideration, the problems involved 
in quantitatively calculating a diffusion coefficient for a 
specific solute-solvent system are quite complex. A con- 
venient first step towards solving these problems is the 
consideration of the rate at which the solvent atoms them- 
selves move around in the solvent lattice, or the self-diffu- 
sion coefficient of the solvent material. After a solution 
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has been found to the problem of self -diffusion, the effects 
of a solute impurity are treated as a superposition problem. 
The relationship between impurity diffusion coefficient and 
the self -diffusion coefficient is discussed in section C-2. 
C-1 Self-Diffusion 

Introduction : In a pure metal, diffusion operates 

by the vacancy mechanism and the self-dif fusion coefficient 
is determined by the frequency with which an atom will jump 
into a vacant neighboring site, to, and by the probability 
that a given neighboring site is vacant, P^. These experi- 
ments are usually performed by plating a very thin film of 
radioactive tracer on a pure metal, annealing, sectionina 
and using the thin film solution to determine D^ . This 
is called the self-diffusion coefficient of the solute 
in the given solvent, and the experiments are done ideally 
at "infinite dilution." For the infinitely dilute alloy, the 
problem is then to estimate whether and by how much w and 
for a solute atom differs from co and P^ for a solvent 
atom. 

C-1.1 Mathematical Model : A general approach to the 

derivation of the diffusion coefficient is to consider the 
problem as a whole sequence of jumps that result in paths of 
atoms migrating through the lattice rather than merely 
jumps between two planes. This concept is referred to as 
the "random walk" approach. 
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Consider that successive atom jumps are vectors 
r, , r^, etc. If r is the total number of completed jumps 
which an atom makes per unit time, then after a time t, a 
total of n = tr completed jumps will have occurred, and the 
atom will have moved an average distance, R( t) 


R(tJ = (6) 

from its initial position. It may be shown that the diffusion 
coefficient, D, is obtained by 


where 

Now 


D = R(t)^/6t 
FTty2 = |RCt) • R(t) 


(7) 


( 8 ) 


For the case of the fee, bcc, and hep crystal systems, all 
the jump vectors will be of equal magnitude, r. Then 


equation (8) may be written as 


RXtJ = nr + 2r z(n-j) cos e 


(9) 


where cos 6j is the average value of the cosine of the angle 
between the i^^ and the ( i + j atom jump. 

If the direction of each jump is independent of 


all earlier jumps, then cos = 0 and substituting equation 
(9) into equation (7j we obtain 

D = ( 1/6 j rr^ . (10) 

The factor ( 1/6 ) come from the fact that only 1/3 of the 
randomly directed vibrations are in the required direction, 
and then substituted into equation (5) . 
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As an example, consider an fee lattice. For this case r = 
12o)N, where 12 is the number of nearest neighbors, u is the 
jump frequency, and is the fraction of vacant sites. 

With a lattice constant a^ , for vacancy diffusion 

D = , 

while for interstitial diffusion 


D = <5a^w 
o 


14 


with 6 being a geometric constant 

The problem of determining D has now been reduced 
to finding a way of calculating and w . 

C-1.2 Correlation Effects : The assumption of uncorrela- 

ted motion made in the derivation of equation (9) is valid 
for interstitial diffusion in dilute alloys. However, it is 
not valid in the case of vacancy or interst it ialcy diffusion. 
Consider vacancy diffusion, after an exchange between an 
impurity atom and a vacancy all of its neighbors are not 
identical; one of them consists of a vacancy, it is most 
probable that the impurity atom will jump back to the position 

occupied by the vacancy. In other words, the mean square dis- 

—2 

placement for the impurity atom after n jumps, R^, will be 

less than that for a vacancy which took the same number of 
—2 2 

jumps, = nr . The ratio of these two quantities defines 


the correlation factor, f, as 

f = lim IF?/R^ 
n-*-® 1 V 


(U) 
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The correlation factor has been calculated by 
Coinpaan and Haven. They demonstrated that the correlation 
factor could be written as 

f = 1 + cos eT / 1 * cos ^ (12) 

where cos eT refers to the average of the cosine of the 
angle between consecutive jumps. For vacancy diffusion in 
a fee lattiee, where the terms have been previously defined, 

D = a^fujN 
o v 

C-1.3 Calculation of N^: According to classical ther- 

modynamics, a given system tends to assume a configuration 
such that the free energy of the system is a minimum. This 
free energy G is related to the enthalpy, H, temperature in 
degrees Kelvin, T, and entropy, S, of the system by 

G = H - TS (13) 

The entropy of any system may be divided into two 
parts, the thermal entropy and the configurational 
entropy S^. The thermal entropy is determined by the 
number of different ways in which the thermal energy of 
the crystal may be distributed over the possible vibrational 
modes of the individual atoms. is related to by the 

Boltzmann relation, where k is Boltzmann’s constant, 

= k In . ( 14 ) 

The configurational entropy of a crystal is determined by 
the number of different ways in which the atoms of the 
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crystal may be arranged over the available number of lattice 


sites. For the case of vacancies in a crystal, let N be 

a 

the number of atoms in the crystal and n^ be the number of 
vacancies, and let all of the lattice sites be equivalent. 
It can be then shown that this case^^ 


W = 
c 


(N^: n^:) 


C15J 


Again the relation between and is given by the 
Boltzmann relation 

= k 1.1 = k In . 

3. V 


(163 


Now assume the energy to create a vacancy in a crystal to be 

h^. The enthalpy of a single crystal containing n^ vacancies 
is increased by a factor riyh^ over that of a perfect crystal. 
For reasons given below, the thermal entropy will increase 
by a factor n^As^ when n^ vacancies are added to the crystal. 
Thus, from these considerations and from equations (13) and 
(16), the free energy of an imperfect crystal as a function 
of n^ and T becomes, 

G(n,T) = Gpeyfect^*^^ ^ ^v^ ‘ V^t ’ ^ (I TTn ^Tr 

(17) 

In order to find the equilibrium value of n^, we note that 
at equilibrium 9G/an^ = 0. Using this fact, and the approx- 
imation that In xl - X In X for x < 1 , we have from (17)^7 


Note: Lower case letters indicate a property of a single 
crystal, while a capital letter indicates a property of an 
entire crystalline structure. 
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( 18 ) 


n -H f r ^ ~T- 

V a eq a eq 

Here, C^y/Na)gq equilibrium atomic fraction of 

vacancies in the crystal. In terms of the equilibrium mole 
fraction of vacancies N^, equation (18) may be written 

ASt -H„ 

^v " ~r“ “ivf • 

In equation (19) , and refer to the thermal entropy and 

enthalpy changes of the entire crystal due to the addition of 
vacancies. 

It should be noted that the foregoing derivation 
assumed that s. and h, are independent of n /N , which 

tv V 3 

would be true if n /N„ , or N,,, were very small. Experiment- 

V a V 

“4 18 

ally, has been fouind to be about 10 in pure metals, 
so this assumption is a reasonable one. 

The enthalpy of a crystal should change with the 
addition of a vacancy, and this energy change equals approxi- 
mately the energy necessary to remove an atom from the inter- 
ior of the crystal and place it on the crystal surface. The 
reason for the change in S^, however, is more subtle. One 
way to estimate this change is to assume an Einstein model 

for the solid consisting of N atoms equally distributed over 

19 

3N harmonic oscillators. When v is the Einstein frequency 
of the atoms in a perfect crystal and hv << k, where h is 
Plank’ sconstant , and k is Boltzmann's constant, then S^, the 
thermal entropy can be shown to be 
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( 20 ) 


S 


t 


3N,k 

a 


[1 + In 


kT| 


for that perfect crystal. In a crystal containing n^ 

the atoms around a vacancy will have a vibra- 
tional frequency v' less than v because the atomic restoring 
forces are reduced- Again using the Einstein model, each 
atom neighboring a vacancy is assumed equivalent to 3 
harmonic oscillators of frequency v’ . For atoms and 
vacancies in a crystal, we have 

3n^x oscillators of frequency v 
(3N^ - 3n^x) oscillators of frequency v' 
where x is the number of nearest neighbor atoms around each 
vacancy. Applying this to equation (20) , we have that the 
thermal entropy of an imperfect crystal is: 


S. = 3n xk |1 + In ^1 + (3N^ - 3n x)k |1 + In | 
t V ' hv ' ^ a V ‘ hv ' 

( 21 ) 

The change in thermal entropy due to the addition of n^ 
vacancies to a perfect crystal can be estimated by subtract- 
ing (20) from (21) and dividing by n^ to obtain 

= 3xk in (v/v') . (22) 

Thus the change in thermal entropy is due primarily to the 
change in vibrational frequency of the atoms around the 
vacancies . 
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Calculation of uj : While the calculation of is a 

straightforward one, accurate determination of the atomic 
jump frequency for a tracer atom is quite difficult. In 
fact, the accuracy of the present calculations of w is so 
poor that such calculations provide no real verification of 
experimentally determined diffusion coefficients. However, 
the following relatively simple theoretical treatment of 
0 ) yields a numerical value of this parameter within at 
least an order of magnitude, and it further provides some 
insight into the temperature dependence of this quantity. 

First, assume a simplified sequence of atomic motions 
required for a tracer atom to jump from one lattice site 
to another, as shown in sequence Figures 2a, 2b, and 2c. 
Before such a sequence can occur, the tracer atom must 
attempt a jump in the direction of the vacant site, with 
the two restraining atoms simultaneously moving far enough 
apart to let the tracer atom slide between them. The 
atomic configuration shown in Figure 2b is called an 
activated complex , and consists of a region containing a 
diffusing atom between two equilibrium sites. The number 
of atoms diffusing per second in one mole of solvent 
material equals the mole fraction of activated complexes 

multiplied by the average velocity of the diffusion atoms 
moving between the restraining atoms v, divided by the 
width of the barrier created by these restraining atoms 6. 
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Free Energy Variation of Lattice During 
Sequence of Atomic Motions Shown in a, b, and c 


MIXi^ANTSM FOR AN ATOMIC JUMP INTO A VACANCY 

Figure 2 
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The average jump frequency per atom, w , becomes 


u = N v/ 6 (23) 

m 

To find an approximate value of N^, it is assumed 
that diffusing atoms enter, remain at, and leave the acti- 
vated complex slowly so that the remainder of the lattice 

. . 21 

has time to continually relax to an equilibrium state . 

Using this assumption, the work done during the resulting 
reversible, isothermal process of moving an impurity atom 
from one site to another is just equal "to the change in the 
Gibbs free energy for the region, G^. From these assumptions, 
at any temperature, there will always be an equilibrium 
number of activated complexes. From the definition of G^, 
in equation (13) 


G = H - TS 
mm m 


(24) 


Since the lattice was assumed to be in thermal equilibrium, 

G will have all of the properties of G,, of equation (17) . 

By a process identical to the derivation of equation (19) 

22 

from equation (17) , we can obtain the equilibrium mole 
fraction of activated complexes at any instant: 




= exp 
m 


m 


TaS. 


m 


-aG 


RT 


exp 


m 


RT 


(25) 


In equation (23), a dimensional analysis will reveal that 
v/ 5 is a frequency. This frequency is conveniently taken 
to be the mean vibrational frequency of an atom about its 
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equilibrium site, v, and represents the frequency with which 
an atom attempts to change sites. From this and from 
equations (23) and (25) we find 


oj = V exp 



+ TAS^ 
m 


RT 


(26) 


It has been determined experimentally that the 
activated complexes do not exist long enough for much of the 
surrounding lattice to reach equilibrium. Thus, it would 
seem that the basis for the derivation of equation (26) is 
open to serious question. However, it does yield an order 
of magnitude agreement with experiment. 

C- 2 Relationship Between Self-Diffusion and Dilute 

Impurity Diffusion 

The problems involved with impurity diffusion are 
more complex than those of self-diffusion, since the atoms 
of both solvent and solute migrate simultaneously. The mi- 
gration of solute atoms is influenced by that of solvent atoms, 
and the presence of impurity atoms in the lattice, changes 
the energy conditions, and thus influences the rate of self- 
diffusion . 

C-2.1 Correlation Effects : In dilute alloys, the 

vacancies are partially bound to the solute atoms, so that 
there is in addition to the correlation between the jumps of 
the solute atoms a correlation between the successive jumps 
of the vacancies. Thus, it is no longer possible for to 
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equal N^, although these two quantities are still closely- 
related. Also, the correlation coefficient becomes much 
less than unity. In this discussion of impurity diffusion, 
it is assumed that the solute concentration is so small 
that the diffusion coefficient of the solute in the solvent 
lattice, or the impurity diffusion coefficient, is inde- 
pendent of the solute concentration. In analyzing the 
behavior of a solute atom-vacancy pair in a solvent lattice 
under these conditions, the frequency with which the solute 
atom and vacancy exchange positions may be taken as the 

frequency with which the vacancy exchanges places with a 
nearest neighbor atom to the solute atom as 0J2, and the 
frequency with which the vacancy exchange places with a 
solvent atom which is not a nearest neighbor to the solute 
atom as Thus is proportional to the frequency at 

which the solute atom- vacancy pair dissociates. Now if the 
vacancy is tightly bound to the solute atom, will be 
nearly zero, and will be much greater than In other 

words, if a vacancy exchanges places with a solute atom 
with a given jump, the probability that it will reverse 
that exchange with its next jump is nearly unity. If this 
is so, the vacancy movement is nearly completely correlated, 
and the solute atom will move through the lattice only as 
fast as the vacancy exchanges with the solvent atoms. Thus, 
the impurity diffusion coefficient for small solute con- 
centrations, Dj( 0 ), will not be given by 
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but by 


DjCO] 'v. 


Dj(0) ^^^o“l’’v 


(27) 


The relationship between and is discussed on page 35. 

In the above equations f<<l, and the correlation effect be- 
comes quite important. Note, the relation between the cor- 
relation coefficient, f, and the quantities , UI 2 , and is 
a function of the lattice geometry. As an example, we can 
consider a face centered cubic lattice, where the expression 
for f is given by^"^ 

(D- + 7o),/2 

^ j I 


f == 


+ oT^ + 7 ( 1 ) ^7^2 


(28) 


Here the return of the dissociated vacancy is ignored. The 
expression for Dj(0) thus becomes 


a 0 ), (uj^ + 7o),/2) 
= i : 7.J2 


(29) 


Present day calculations of ^2 and co^ yield ambiguous results 

so the best practical approach is to estimate the value of 

25 

Dj(0) for several limiting cases. These are: 

(1) the rate of dissociation of solute atom vacancy pairs 

is much less than the exchange rate between the vacancy 
and the nearest neighbor to the solute atom, or 


'3' 


Then we have from from equation (29) 


Dj(0) '^.a; 


2,“2“l*’v 


“l “2 


(30) 
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( 2 ) the vacancy-solute atom exchange rate is much more 

rapid than the two vacancy-solvent atom exchange rates, 
or the vacancy-solute atom pair is tightly bound. Then 

w^>>(w 2 + in (28), and we have 

DjCO) a^(u)2 + 70)3/2)?^ ( 31 ) 

(3) the solute atom jumps slowly, or o)j^<<(tjj 2 70 ) 3 / 2 ) in 

equation (29'). Thus 

DjCO] •V (32) 

(4) the "solute" is a solvent tracer, so that = 0)2 " <^^3 = 
0 ). Then, since here = N^, and D^(0) = Dj(0) 

D^CO) = 2£a^N^io . (33) 

It must be emphasized that the above results are 
for a face centered cubic lattice. Other lattice geometries 
will give results of essentially the same form, but with 
different numerical constants. 

C-2.2 Valence Effects : Due to the electronic effect of 

an impurity atom with a valence different than that of the 
surrounding host lattice a vacancy is often attracted to 
the solute atom. Consider the case of a divalent impurity 
in a monovalent lattice. According to the free electron 
theory of metals the free electron density in the immediate 
environment of the impurity will be greater than that 
throughout the lattice as a whole, but it will not be great 
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enough to completely neutralize the excess positive charge 
of the solute ion. This ion is said to be "screened,” and 
it is surrounded by an electrostatic potential of the form 

V(r3 = (Ze/r ) exp (-qV ) . (34) 

Where r is the radial distance from the impurity ion, Z is 
the number of excess electrons per impurity ion and q is 
the screening parameter, which gives a measure of the 
neutralization of the excess charge of the impurity ion. 

Let a solvent atom be removed from a site adjacent to the 
solute ion, and make the assumption that the electron dis- 
tribution throughout the lattice will not be appreciably 
distributed by this removal. If this is the case, the 
vacancy resulting from this removal will have an effective 
charge of -e. Its energy on the site next to the impurity 
ion then will be reduced by 

eV(r) - E(r) = ( Ze^/r ) exp ( -qr ) (35) 

as compared to a site surrounded by normal solvent atoms. 
Since the energy of the vacancy is lower next to the 
impurity ion, the concentration of vacancies on any nearest 
neighbor site will be increased to 

Py “ ^v )/kT) . (36) 

This effect will also decrease on a site next to a solute 

atom by E(r ) as compared to on a site surrounded by 
solvent atoms. This decrease in potential energy also makes 
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for the solute less than for the solvent. Finally, 
there is a force e(dV/dr) tending to draw the vacancy and 
the impurity ion together. If the impurity ion and the 
vacancy begin to move together, this force tends to pull 
them together until is decreased by roughly r^e (dV/dr) 

( 

C-2.3 Size Effects : The effects on the impurity dif- 

fusion coefficient of the size of the solute atom relative 

2 7 

to that of the solvent atom have been considered by Swalin 
28 

and Overhauser. 

Swalin assumes that the solute atoms are com- 
pressible spheres and the solvent lattice is an elastic con- 
tinuum. For example, if the solute is larger than the 
solvent atom, will be increased by the additional strain 
required at the saddle point. This change in is a 
balance between the strain on the diffusing solute atom and 
the lattice continuum. Although it would seem that a large 
solute atom would elastically strain the lattice and thereby 
attract vacancies and lower H^, Swalin concluded that the 

actual change in is negligible to a first order approxi- 
27 

matron. 

28 

Overhauser considered the effects on the 
impurity diffusion coefficient of the presence of inter- 
stitials in the solvent lattice. Since these interstitials 
dilate the lattice locally, the average interatomic dis- 
tance between the matrix atoms will be increased. This 
lattice expansion will reduce the work required to squeeze 
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an atom through the saddle point, so should be decreased 
by the addition of an interstitial solute whose atoms are 
larger than those of the solvent. Likewise, the addition 
of an interstitial solute which contracts the solvent 
lattice should increase 

Although the impurity diffusion concept has been 
studied in detail for many elements, no single correlation 
between diffusion parameters D^, the frequency factor, and Q, 
the activation energy and the properties such as valence and 
size has been found to fit more than one or two elements. 
C-2.4 Solute Concentration Effects : When the solute 

concentration becomes greater than some threshold value, 
Cwhich varies for different systems) , the solute atoms begin 
to interact and the diffusion coefficient for the solute as 
well as the self-diffusion coefficient for the solvent 
material becomes a function of the solute concentration. 

In order to obtain a qualitative concept of the nature of 
these functions, we consider the previously presented 
equations for the self-dif fusion of the pure solvent^ C33) 
and the dilute impurity diffusion coefficient for the solute, 
in a face centered cubic solvent lattice, (28), i,e., 


and 


Ds( 0 ) = 2fa;.^N^ 


% 1^2 

DjCO) = 2| ^° ^ IP. 


1 " "2 " 


TJITTTz 


(33) 

(28) 
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Now i£ Py = N^ and ^.102, - then no matter how large 

is, D^( 0 ) will about equal Dj( 0 ). However, i£ Dj( 0 ) is 
greater that DgCO) very dilute alloys, it £ollows that 
the impurity atoms attract vacancies so that > N^. Then 
the solvent jump £requencies ^2 and are increased near 
a solute atom. But if P, > N . or near a solute 

V V 4 0 

atom, then is also greater near a solute atom. As the 

solute concentration increases, so does the number of affected 

solvent atoms. Thus, if Dj( 0 ) < D^( 0 ), an addition of solute 

will increase . Similarly, if Dj(Oj > D^( 0 ), an addition 

of solute will decrease . 

s 

To provide an estimate of the magnitude of the 
change in the self-diffusion coefficient, , with solute 

concentration, the following empirical equations have been 
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derived for solute concentration _< 0.01: 

= ( 1 - ANj )( D^(0)) + A3 NjDj (37) 

Dj = Dj(0) (38) 

Here, Nj is the solute concentration, A is the number of 
nearest neighbor sites around each atom, and 3 is an arbi- 
trary constant ^ 1.0 which must be experimentally determined. 
Thus, Dj does not vary much with solute concentration for 
solute concentrations less than 0 . 01 . For solute concentra- 
tions greater than 0.01, the relation between the impurity 
diffusion coefficient and solute concentration is given by 

Dj = Dj(0)( 1 + yNj ) (39) 
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Here, y is a constant ^1.0. The values of the constants 3 
and y vary from system to system. Specific theoretical 
calculations of their values are quite difficult and are 
not usually attempted. 

D. Summary of Equations for Diffusion Coefficient 

The following comprise a summary of the most use- 
ful previously derived diffusion relationships which are 
most useful in comparing experimentally observable parameters 
with theoretical models for the diffusion process. 

For vacancy self -diffusion in a pure solvent with 
correlated tracer jumps, the self -diffusion coefficient, , 
is given by 

D^(0) - 2fa^N^co (40) 

where f is the correlation coefficient, is the equili- 
brium vacancy concentration, oj is the tracer jump frequency, 
and a^ is the lattice constant. For dilute impurity vacancy 
diffusion, the diffusion coefficient of the solute in the 
solvent lattice is 

Dj(0) = (41) 

Here is the frequency with which a solute atom changes 
places with a nearest neighbor vacancy , and is the 
probability that a nearest neighbor site to the impurity 
atom is vacant. Because of attractive forces between the 
vacancies and solute atoms have already been discussed, 
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An estimate o£ the 


in equation (27) will not equal N^. 
value o£ in terms of for a given solvent-solute sys- 
tem may be obtained from equations (33), (34), and (35). 

The coefficient of correlation £ is a function of lattice 
geometry. For self-diffusion, f lies between 0.50 and 0.79;^^ 
for a face -centered cubic solvent structure, £ for impurity 
vacancy diffusion is given by 


f = 


^2 ^ 7 ( 0^/2 

<iJi + a>2 7 u)2 


(42) 


Here is the frequency with which a vacancy exchanges 
position with a nearest neighbor atom to a solute atom, and 
^2 is the frequency with which a vacancy exchanges position 
with a solvent atom, which is not a nearest neighbor to a 
solute atom. 

From equilibrium thermodynamics, the equilibrium 
vacancy concentration in a crystal was found to be 


AS. 


H. 


Ny exp 

In a similar manner, an equilibrium approximation allowed 


(43) 


us to arrive at an expression for 

+ TAS^ 
m m 

= V exp 


(44) 


Combining equations (19), (26), and (40) allows us to write 

for the vacancy self-diffusion coefficient in a pure or 
very dilutely contaminated solvent 
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- ASj^ + AS ”^f ” ^ 

Dg(0) = 2£a^ V exp | ^ |exp | rt^I * 

Likewise, combining equations ( 19 ), ( 26 ) , (27) , and (35) , we 
find for dilute impurity vacancy diffusion 

AS. + AS^ -H. E(r ) 

Dj(0) = 2fa^ V exp| ^ exp | — ^ — -\ exp| — j^l • 

(46) 

Here E(r ) may be estimated from equation (34), and f for a 
fee lattice is given by equation (28). Generally, f will be 
given in terms of the various jump frequencies described 
above. Of these frequencies, only can be calculated with 
reasonable certainty; because of uncertainties in the values 
of 0^2 and the usual practice is to consider the behavior 
of Dj^(O) in the limiting cases given in equations (30)^ (31), 

and (32). 

Equations (37), (38), and (39) give the variation 
in the solvent self-diffusion coefficient D^(0), and vacancy 
impurity diffusion coefficient, Dj(0), with solute concentra- 
tion. Finally, for interstitial diffusion, the diffusion 
coefficient, D, is given by 

D = a^ V exp |-^| exp |~^| . (47) 

Thus, interstitial diffusion can be considered a much simpler 
process than vacancy diffusion. For the present, investigation 
it is probable that both vacancy and interstitial diffusion 
mechanisms take place. 
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Equations C45) , C46) , and (47) are useful in cal- 
culating values of D from theoretical models of the diffusion 
process. To check the validity of a given model, a compari- 
son is made of the calculated values of the factors compri- 
sing D with corresponding experimentally determined values. 
Several workers have devised experimental methods for deter- 
mining the individual terms of equations (45) , (46) , and 
(47), such as aSj., AS , H , etc.. However, the most common 
procedure is first to express these equations in the form 

D = exp ( -Q/RT ) . (48) 

Then, after determining the diffusion coefficients of the 
system under consideration at different temperatures, Q may 
be found from the slope of the plot of the natural logarithm 
of D versus ( temperature , °K ) while the intercept of this 
plot with the ordinate axis yields D^. The quantity Q is 
called the activation energy, and is the sum of the enthalpy 
terms in each of the equations for D, (45), (46), and (47). 
This quantity is thus associated with the energy necessary 
for the solute ion to make a single jump in the solvent lat- 
tice. The symbol, D^, called the frequency factor, contains 
the correlation coefficient, jump attempt frequency, lattice 
parameter, and the entropy changes associated with an atomic 
jump. Hence, gives a measure of the disorder induced in 
the solvent lattice by an atomic jump. 


35 



The " random walk '* is not the only approach for 

the explanation of the phenomena of diffusion. Some of the 

more interesting theories are presented by: Feit discusses a 

dynamical theory, Huntington and Seitz with a discussion 

32 

from the modern theory of metals, Lazarus discusses the 

33 

effect of ’* screened " electrons on diffusion, Johnson 

34 

presents a " hole ” theory for diffusion, while Zener 

35 

presents a model for ring diffusion. All of these studies 

present arguments which are valid for some very limited 

cases and were therefore not used in this work. For further 

36 

reading on diffusion, the studies by LeClaire , Wert and 

37 38 39 

Zener, Turnbull and Hoffman, and Zener are suggested. 
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Part III 


MATERIALS AND APPARATUS 

This research depended upon mass spectrometr ic 
instrumentation o£ high sensitivity. It was also condition- 
ed by the type and form of samples which were amenable to 
diffusion analysis by mass spectrometr ic means. 

The samples of zircaloy-2, tantalum, and 304 
stainless steel ribbon were obtained from H. Cross Company 
of Weehawken, New Jersey, with a detailed analysis of the 
composition of each sample material. These appear in Tables 

I, II, and III. 

The mass spectrometer used in this research was 
a 90^ double focusing instrument with a resolution of 
approximately 1/1000, designed and built by F. A. White and 

J. C. Sheffield. 

A. Description of the Filament 

A single filament thermal ionization source. 

Figure 3, was used in all experiments. Each of the three 
sample materials (304 stainless steel, zircaloy-2, and 
tantalum) were used as the filament material. The ribbons 
of the zircaloy-2 and tantalum were 0.030" (0.076 cm) x 
0.001" (0.0025 cm), while the 304 stainless steel ribbon 
was 0.030" (0.076 cm) x 0.0005" (0.0013 cm). 

Heremetically sealed, insulated feedthrus, on 
which the filament ribbons were spot-welded, were obtained 
from Electrical Industries, Murray Hill, New Jersey. The 
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Table I 

COMPOSITION OF 304 STAINLESS STEEL 


Carbon 

. 056 

Manganese 

1,65 

Phosphourus 

.018 

Sulphur 

.012 

Silicon 

.49 

Chromium 

18.53 

Molybdenum 

. 31 

Copper 

.25 

Nickel 

9.14 

Cobalt 

. 03 
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Table II 


COMPOSITION OF TANTALUM 


TANTALUM 

Atomic Number 

73 


Atomic Weight 

180.95 

Ta 

Crystal Structure 

BCC 

Melting Point ®C 

2996 

Density gms/cc (20®C) 

16.6 


Ibs./in^ 

0.600 


MARZ GRADE 


MATERIALS ANALYSIS 


Analysis: ■ Mass Spectrographic ■ Vacuum Fusion (Gases) □ Emission Spectrographic ■ Conductometric (Carbon) 


Element 

Content (ppm) 

Element 

Content (ppm) 

Element 

Content (ppm) 

Element 

Content (ppm) 

0 

< 5.0 

S 

2.0 

Cu 

0.3 

Sn 

0.2 

c 

< 5.0 

Cl 

2.0 

Zn 

< 0.1 

Te 

< 0.2 

H 

< 5.0 

K 

0.1 

Ge 

< 0.4 

1 

< 1.0 

N 

< 3.0 

Ca 

0.4 

As 

< 0.1 

W 

3.0 

B 

< 0.1 

Sc 

< 0.1 

Br 

< 1.0 

Pt 

< 5.0 

F 

< 1,0 

Ti 

< 1.0 

Zr 

< 1.0 

Au 

< 5.0 

Na 

< 0.4 

V 

0.2 

Nb 

20, n 

Pb 

< 0.6 

Mg 

< 0.2 

Cr 

3.0 

Mo 

< 0.5 

Bf 

< 0.3 

Al 

0.5 

Mn 

< 1.0 

Ag 

< 0.4 

Co 

< 0.3 

Si 

0.7 

Fe 

iOO 

Cd 

< 0.3 

All 


P 

0.3 

Nl 

1.0 

In 

< 0.2 

others 

< 0.1 


METHOD OF PREPARATION: Electron Beam NOMINAL PURITY (TOTAL): 99 996% 

Zone Refined 

STANDARD SINGLE CRYSTALS 


$/in. (1" Minimum Order) 

diameter H" (.450") 

$100. S125. $150. 5225. 


Random Orientation. Oriented crystals, $75. additional orientation charge per crystal, independent of length. (100), 
(110), and (111) available. 

SPECIAL APPUCATION CRYSTALS 

Cut from */ 2 " nominal diameter (elliptical cross-section) to approximately Vs" thickness, one face polished strain-free 
to expose specified crystallographic plane: 

Within 1° $6W./slice; additional slices of same orientation $200. each. 

Within 4° $325. /slice; additional slices of same orientation $150. each. 

POLYCRYSTALUNE RODS 


< 100 gms 100-300 gms 1-3 lbs. 

diameter in./lb. gms/in. $/gm $/gm S/lb. >3 lbs. 


Ks" 

137 

3.32 

$3.00 

$1.70 

$550. 

$500. 

y," 

34 

13.27 

2.00 

1.60 

500. 

450. 

X" 

15 

29.84 

1.85 

1.55 

475. 

425. 

y 

8.5 

53 

1.75 

1.45 

450. 

400. 

1" 

2.1 

2.20 

- 

1.38 

400. 

350. 


>lb. on request. 
Minimum Order; $ 50.00 


39 



TABLE III 


COMPOSITION OF ZIRCALOY-2 


ZIRCALOY-2 INGOT 
HEAT NO. 387266 


COMPOSITION IN PERCENT 



Top 

Bottom 

Sn 

1.45 

1.48 

Fe 

0.15 

0.16 

Cr 

0.11 

0.11 

N1 

0.05 

0.06 

Zr 

B A L A N 

C E 


IMPURITIES 

IN PPM 

A1 

47 

49 

B 

0.2 

0.2 

C 

110 

100 

Cd 

<0.2 

<0.2 

Co 

<10 

<10 

Cu 

16 

15 

H 

7 

<5 

Hf 

41 

39 

0 

1180 

1150 

Mn 

<25 

<25 

N 

46 

38 

Si 

73 

79 

Ti 

<25 

<25 

W 

<25 

<25 

U 

1.5 

1.8 

V 

<25 

<25 


INGOT HARDNESS 

, BHN 


Range 170 - 179 

Average 176 
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SINGLE FILAMENT THERMAL IONIZATION SOURCE 
FIGURE 3 
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feedthrus are designated as El type A-40W-SS Mod H. 

B. The Mass Spectrometer 

The mass spectrometer separates chemical elements 
according to their mass to charge ratio. The mass spectro- 
meter shovm in Figure 4 was used exclusively in this study. 

A schematic of this two-stage, double focusing instrument 
is shown in Figure 5. The major components of the instru- 
ment include: 

(a) Ion source 

(b) Electrostatic analyzer 

(c) Magnetic analyzer 

(d) Beam defining slit system 

(e) Ion detector and counting system 

(f) Vacuum system 

Thermal Ionization . For mass spectrometric anal- 
ysis, some method is needed to vaporize and ionize the sam- 
ple into a beam of positively charged ions before entering 
the electrostatic and magnetic filters. For the analysis of 
metals, thermal ionization is an effective means of produc- 
ing gaseous ions especially if the metal has a low ioniza- 
tion potential compared to the filament’s work function. As 
the sample vaporizes, some atoms or molecules will leave the 
surface of the filament with the loss of one or more of their 
electrons. The hot filament surface has a higher affinity 
for electrons than the metal itself. This is due in part, 
to the intrinsic nature of the metallic crystal along with 
the low binding energy of the outer shell electrons of most 
metallic atoms. 
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TWO-STAGE MASS SPECTROMETER 


FIGURE 4 



amplifier /DISCRIMINATOR / DRIVER 



Two Stage, Double Focusing Mass Spectrometer 


Figure 5 
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The probability o£ an atom evaporating from a sur- 
face in the ionic state rather than as a neutral atom is 
expressed by the relation. 


n_ . exp . exp (49) 

n 

where n'*’ is the number of positive ions, n° is the number of 
neutral specie, e is the electron charge, IP is the ioniza- 
tion potential of the sample in electron volts, cj) is the 


surface work function of the filament material in electron 


volts, k is Boltzmann's constant, and T is the surface 
temperature in °K. 

It can be seen that for elements where IPxfi the 
relative number of ions to neutrals produced will increase 
with temperature. 


Ions resulting from thermal ionization on the 
filament surface are accelerated to high energy and defined 
by collimating slits in the ion source. These high energy 
ions pass through the electrostatic analyzer where they are 
energy filtered. Ions with a discrete kinetic energy 
Caccelerat ing voltage) including a small energy spread 
pass through the center defining aperture. The defining 
slit determines the energy spread accepted. This filtered 
beam of ions then enters the magnetic analyzer where the 
beam is momentum filtered. Those ions of appropriate 
energy and proper momentum (essentially mass) come to a 
focus at the final slit and strike an electron multiplier. 
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Current pulses result from ion-electron conversion and sub- 
sequent secondary electron amplification in the electron 
multiplier tube. Output pulses are further amplified and 
counted by an amplif ier/discriminator/scalar system. Thus, 
ions of the same mass are individually counted and the 
exceedingly small ion currents can be quantitatively 
measured . 

By properly changing the accelerating potential 
of the ion source, ions of other mass numbers are focused 
on the electron multiplier. By comparing counting rates 
for each mass position the isotopic ratio can be quantita- 
tively determined. 
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Part IV 


EXPERIMENTAL METHOD 

A. General Experimental Techniques 

From the literature, the most common method for ob 
taining diffusion coefficients is that of serial sectioning. 
However, a mass spectrometric approach was chosen for this 
work. The basic concept of the mass spectrometric technique 
for determining diffusion coefficients is that a correlation 
exists between the solute atoms moving through a solvent 
lattice and the fraction of these atoms which are evaporated 
as ions, from the crystal surface at a given temperature. 

In this work, a thin ribbon filament of solvent 
material is used with its own impurity or solute concentra- 
tion. No doping of the solvent was required and thus there 
was no alteration of the crystal lattice. Since the mass 
spectrometer source chamber pressure is well below the vapor 
pressure of the solute material in the temperature range con 
sidered, the solute atoms diffuse to the filament surface, 
and evaporate immediately. This results in a zero solute 
concentration at the filament surface after the filament 
temperature is raised to the desired level. Further, as the 
ribbon thickness is much less than its length or width, the 
diffusion problem may be considered to be one dimensional. 

Application of Fick’s Second Law for Diffusion 
in a Thin Plate : Since the solute concentrations 
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under consideration are very small, the diffusion process 
occurring in the ribbon may be described by 

J(x,t) = -D (503 

and 

3 C Cx , t ) _ _ j j3 C (x , 1 3 ( 5 1) 

3t 3X^ 

where J(x,t) is the solute flux, C(x,t) is the solute con- 
centration, and D is the diffusion coefficient. For the 
experimental conditions previously described, the solution 
to Pick's second law for impurity evaporation from a thin 
plate yields the solute concentration, as a function of 
position in the filament and evaporation time,t. 

Solving equation (51) by the method of separation 

of variables, the solute concentration may be found to be 

2 

C(x,t) = C A sin Ax + B cos Ax ) exp (-A Dt) (52) 
where = - (1/dT) (dT/dt) or - (1/X) (d^X/dx^) , D is the di£- 
fusion coefficient, t is time, and A and B are constants of 
integration. 

For solute depletion from a thin solvent plate with 
an initial solute concentration f(x) across its thickness, 
the boundary and initial conditions become 

C - 0 at t ^ 0, for x = h and x = 0 
C = f(x) at t = 0, for h > x > 0 
With these conditions, equation (52) becomes 

C(x,t) = ^ E sin ^ exp / f(x') sin dx' 
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Where h is the filament thickness. Thus from Pick’s first 
law, equation (50), the solute atom flux at the filament 
surface ( x =0 ) becomes 

n 2 2 

J( 0 ,t) = ( 2 /h)Z exp I ^ I /f(x') sin dx’ (54) 

h^ h 

or more conveniently, where n is an integer, 

2 2 

J(0,t) = A Z expl ^55) 

2 

where A is a constant. For h < 16Dt the error in using 
only the first term is less than 11 .^^ 

Positive Ion Emission from a Hot Metal Filament : 

The combination of high temperature and solvent work function 
causes a constant fraction of the evaporating solute atoms 
to ionize at the filament surface at a given temperature. 
According to the Langmuir- Saha equation, the ratio of the 
number of emitted positive ions (n^) to emitted neutrals (n°) 
is 

n‘^/n° = B exp { ( 4> “ IP)/kT} ( 55 ) 

where B is a proportionality constant, <}> is the filament 
work function, IP is the ionization potential for the solute, 
k is Boltzmann’s constant, and T is the temperature of the 
filament in degrees Kelvin. 

The mass spectrometer accelerates, mass analyzes, 
and detects these emitted solute ions while excluding all 
other ionic species from its detector. By combining 
equations (54) and (55), and assuming that the atomic trans* 
port time is greater than h /16Dt, the time variation of the 
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detected ion current is 

J(0,t) = A exp C -ir^Dt/h^ ) . (57) 

B. Experimental Procedures 

The procedures used in this work are detailed 
below to allow reproduction of the experiments. A general 
procedure followed at all times was the wearing of lintless 
gloves when handling anything going into the vacuiom environ- 
ment. The procedures described are: 

1. Cleaning 

2. Spotwelding 

3. Mass spectrometry 

4. Calibration of mass spectrometer 

5. Filament temperature measurement 

6. Data acquisition 

7. Treatment of data 

8. Error analysis 

1 . Cleaning : 

Filament holders: Clean metal surfaces with 4M HNO,, 

and rinse with de-ionized water. 

Dry under the infrared lamp. 

Filament legs : Place in an alcohol bath immersed 

in the ultrasonic cleaner. Rinse 
with de-ionized water and dry under 
an infrared lamp. 

2 . Spotwelding : 

By visual inspection, make sure electrodes are clean 
Set to 18-25 watt-seconds 

Place one end of filament next to the supporting legs 
and then spotweld ^ 

Do the same for the other end of the filament and leg 
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3. Mass Spectrometer : 

Opening Source Region: 

Close main gate valve, isolation valve, and the 
roughing valve 

Open the foreline valve 

Connect up nitrogen gas, turn on, and let it into 
the source chamber 

Let chamber come to atmosphere, put clips on 0-ring 

Before dismantiling anything, make sure power is OFF 

Screw filament assembly onto the top plate of the 
ion source. 

Closing Source Section and Vacuum Pump Down: 

Make sure no lint is on 0-ring or glass , and then 
put glass back on 

Turn nitrogen gas off and close inlet valve 

Open roughing valve, let it run for a few minutes 

Close roughing and open foreline valve 

Turn on first ionization gauge 

Make sure liquid nitrogen is in cold trap 

Open up the diffusion pump valve, do this SLOWLY 
and watch pressure gauge 

Turn on automatic ionization gauge 

Opening System and Operation: 

Open isolation valve to system, after chamber 
pressure has reached 5 x 10~' torr 

Set electron multiplier tube to -4500 volts 

Close gate valve on multiplier tube’s ion pump 

Turn on polarity switch to minus, on the electro- 
static lens power supply 

Turn on main power to mass bridge C ion source 
voltage divider ) 
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Turn voltage control on filament heater, adjust 
current control for the temperature setting 

Turn up voltage to desired setting, CAREFULLY 

Begin data acquisition 

System Shut Down: 

Turn all voltage settings OFF 

Close isolation valve 

Open ion pump to electron multiplier tube 
Turn multiplier tube voltage to -500 volts 
Close diffusion pump gate valve 
Turn OFF all power 

Note: Liquid nitrogen will evaporate away 
4 . Calibration of the Mass Spectrometer : 

It is necessary to calibrate the mass spectrometer 
before undertaking any experiments measuring diffusion coef- 
ficients. Using equation 

= (144) (58) 

where B is the magnetic field in gauss , R is the radius of 
curvature (in cm.) of the magnetic analyzer, m is the mass 
of the isotope of interest in amu. , and V is the accelerating 

O' 

potential in volts, the appropriate magnetic field-accelerating 
potential values for a given mass can be determined. 

Since the alkali metal atoms were to be examined, 
calibrations were made using filaments which had solutions of 
either rubidium or cesium evaporated on them. Rubidium has 
two naturally occurring isotopes, mass 85 and mass 87. De- 
tection of one mass and calculation of the magnetic field 
allowed the determination of the accelerating potential of 


52 



the other isotope with the same magnetic field. In practice, 
an accurate calibration requires many repetitive measurements. 
Cesium-133 is the only stable isotope of cesium, and provides 
an accurate determination of the magnetic field. 

Ion implanted, as well as micro -pipetted filaments 
of Cs-133 were also used for calibration purposes. The iden- 
tical method was followed. 

5 . Filament Temperature Measurement : 

Using a precision optical pyrometer, filament 
temperatures were determined by visually matching the color 
of the filament with the color of a calibrated, resistively 
heated wire within the pyrometer. It is accurate to 
within ±3°C. 

6 . Data Acquisition Procedure : 

The mass spectrometer data was obtained as follows : 

1. The sample filament was placed in the ion source 
and the source chamber pressure was evacuated 

to a pressure of about 10 torr. 

2. The source voltage was set at the accelerating 
potential of the desired diffusing isotope. 

3. The filament was raised in temperature until 
ions were detected at the given accelerating 
potential. The filament temperature was then 
determined. 

4. The ion emission rate as a function of time was 
monitored until a slowly decreasing emission 
rate was observed. The latter indicated that 
the maximum emission rate had occurred and the 
exponential decay portion predicted by diffusion 
theory, had been reached. 

5. At this point, the ion current was recorded for 
one second every four seconds utilizing a 

frequency-period-counter which has an accuracy 
of ±1 count. 
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6 . 


The procedure was repeated for successively 
higher temperatures . 

7 . Treatment of Data : 

A graph of the natural logarithm of the impurity 
ion current (represented by counts per second) 
from a particular solvent lattice versus time, 
for each temperature was plotted 

The slope of the " linear best fit" line through the 
data was calculated 

The diffusion coefficient was then calculated by 

2 

( slope of the line )♦( thickness of filament/ir ) 

= -D 

4 o 

A graph was made of the In D versus 10 /T ( K) 

Utilizing a PDP 8/1 computer, the linear best fit 
line to the data points was made. The y-intercept 
of this line is the frequency factor, D , and was 
calculated by the computer. The correlation co- 
efficient of the data to the line was also calculated. 

The activation energy, Q, was determined by multi- 
plying the slope of the line from the line in the 
preceding step by the Rydberg gas constant in units 
of kcal/mole/®K. 

8 . Error Analysis : 

Programs were written to provide the linear best 
fit line from a set of data points. The program also gave 
the detree of fitness in the form of a correlation coefficient. 
All error values were determined by calculating the product- 
moment coefficient of correlation and using a 95% confidence 
level. A final program was written to obtain the standard 
deviation for the activation energies. These programs appear 
in Appendix 1 and 2 . 
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Part V 


RESULTS AND DISCUSSION 

Hollowing the procedure previously indicated, graphs 

of the In (alkali metal ion count rate) versus time were 

obtained. Figure 6 is an example of how the ion count rate 

varies with time. These superposed curves illustrate the 

detection of singly charged, mass resolved ions at several 

temperatures. The counting rates are relative, but most 

'14 

experimental measurements were in the region of 10 to 

- 17 5 2 

10 amperes, corresponding to counting rates of 10 to 10 

ions per second. 

Tables IV, V, and VI show the experimental results 
for the diffusion of the alkali metals in 304 stainless 
steel, tantalum, and zircaloy-2. All five naturally occurring 
alkali metals were found in tantalum and in 304 stainless 
steel; however, lithium was not found in sufficient quan- 
tities to determine its diffusion parameters in the zircaloy-2 
sample . 

Activation energy plots for these diffusion systems 
are shown in Figures 7 thru 20. The frequency factor, D^, 
and the activation energy, Q, were calculated from these 
plots and are presented in Tables VII, VIII, and IX. 

The data derived from the activation energy plots 
can be applied to the metallurgy of the system. For example, 
a sharp discontinuity can be attributed to an allotropic 
phase change in the host material. Since pure zirconium has 
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TIME { 5 SECOND INTERVALS ) 


DIFFUSION OF Rb IN TYPE 304 STAINLESS STEEL 
FIGURE 6 
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TABLE IV 


Cn 


DIFFUSION COEFFICIENTS OF THE ALKALI METALS IN TANTALUM 


Ll-7 


T(°K) 

2 

D(cm /sec) 


2 

D(cin /sec) 

I(^ 


-9 


-8 


1080 

4.52-10 1 

1522 

1.06*10 ° 

1110 

1107 

5.03-10’“^ 

1534 

7.23- 10"^ 

1.23- 10"“ 

1145 

1152 

6.86-10“^ 

1.10-10“^ 

1.91-10"*^ 

1570 

1172 

1179 

1575 

8.37-10"^ 

1.31*10"“ 

1.01-10"“ 

1.34*10"“ 

1.13-10"“ 

1199 

1289 

1605 

1235 


q 

1621 

1247 

1381 

3.08-10 Q 

1650 


1439 

4.37-10 „ 

1650 

1267 

1449 

4.77*10-^ 

1669 

1.34*10"“ 

1.43*10"“ 

2.68*10"“ 

3.76-10" 

1300 

1462 

1481 

7.70-10 1 
7.18-10 q 

1675 

1721 

1314 

1342 

1490 

5.58*10“q 

1802 

1355 

1520 

8.92*10‘^ 



1383 


Na- 

•21 


1403 



1422 

1033 

1058 

7.31*10"q 

7.66*10“q 

1311 

1330 

3.27-10"q 

3.30-10"q 

1447 

1107 

7.94*10 q 

1348 

4.69* 10_q 


1110 

8.01-10"^ 

1377 

4.46-10 g 

1050 

1140 

1.32-10"^ 

1464 

4.71-10"g 

1076 

1176 

1.68-10 p 

1.23-10"“ 

1.66-10"“ 

1.63*10" 

1466 

6.39-10"g 

1083 

1190 

1488 

5.97*10"g 

1111 

1230 

1490 

2.99-10"g 

1133 

1230 

1527 

6.86-10"^ 

1.23-10"“ 

1.31-10"“ 

1.09-10"“ 

1155 



1610 


1238 

2.80*10 1 

1634 

1167 

1289 

2.10*10"q 

1661 

1227 

1290 

2.90*10"^ 



1261 





1366 


K-39 


D(cm /sec) 


2.18-10 

3.47-10 

4.18-10 

4.64-10 

7.44-10 


-9 

-9 

-9 

-9 

-9 

-9 


2.02-10 

2.04*10 

2.51*10 

2.55-10 

3.78-10 

2.86-10 

4.01-10 

4.11-10 

4.36-10 


-9 

-9 

-9 

-9 

-9 

-9 

-9 

-9 

-9 


T(“k) 

D(cm^/sec! 

1471 

4.77-10"g 

1495 

6.78-10"q 

1497 

4.52-10"^ 

1499 

5.49-10 g 

1524 

7.23-10"„ 

1529 

6.06- 10"g 

1577 

9.46*10 g 

1595 

9.04* 10"o 

1616 

1.29*10"® 

1664 

1.26*10"“ 

1672 

1.45-10"“ 

1718 

1.45*10"® 

1745 

2.55-10"“ 

2.00*10"“ 

1767 

1828 

2.64-10"“ 


Rb-85 


2.60-10 ^ 
8.00-10‘q 
6 . 10 - 10 “^ 

9.80- 10”^ 
1 . 10 - 10 "^ 
1.60-10 

5 . 20 . 10 *'^ 

5.40-10"q 

6.30-10"q 

7.80- 10“^ 


1372 7.80-10 
1395 1.40-10’ 
1437 1.20*10 
1443 1.60-10 
1538 2.30-10 
1560 2.70-10 
1582 li90-10 
1616 2.80-10 
1667 2,60-10 
1695 2.80-10 
1783 3.10-10 
1789 4.30-10' 



I 


00 


I 


TABLE IV CONCLUDED 


Cs-133 


T&I 

2 

D(cm /sec) 
-9 

T(°K) 

992 

4.46*10 1 

1377 

1040 

5.50*10 q 

1393 

1070 

7.32*10 q 

1422 

1082 

7.96*10 q 

1506 

1100 

9.58*10“^ 

1543 

1122 

1.24*10"^ 

1555 

1124 

1.03*10“^ 

1600 

1133 

1.01*10"^ 

1623 

1181 

1.38*10"“ 

-9 

1631 

1650 

1178 

3.74*10 

1712 

1289 

5.19*10"„ 

1745 

]316 

8.42*10"q 


1326 

5.33*10 


1366 

7.35*10” 




1.09*10 

9.96*10 

1.20*10 

1.54*10 

1.59*10 

2.20*10 

1.58*10 

2.14*10 

1.75*10 

1.61*10 

1.78*10 

2.47*10 


-8 

-9 

-8 

-8 

-8 

-8 

-8 

-8 

-8 

-8 

-8 

-8 



DIFFUSION COEFFICIENTS 


on 

<£3 


Na-23 


T(°K) 

2 

D(cm /sec) 

T(\) 

914 

5.67-10"o 

1135 

933 

5.20*10-^ 

1149 

955 

1.10*10 “ 

1157 

991 

1.78*10“'® 

1181 

1013 

1.88*10“® 

1185 

1060 

9.18*10"q® 

1208 

1229 

1075 

1.02*10“q 

1266 

1076 

1091 

1.01*10"^ 

1.03*10“^ 

1326 

1117 

3.65*10“^ 

1357 

1121 

2.41*10“^ 

1366 



3.21*10 

4.A2-10 

3.43*10 

7.61*10 

7.66*10 

8.11*10 

9.51*10 

1 . 86*10 


-9 

-9 

-9 

-9 

-9 

-9 

-9 

-8 


5.30*10 

1.23*10' 

2.03*10' 


0^ 00 <X) 


V 


THE ALKALI METALS IN ZIRCALOY-2 


K-39 


T(°K) 

2 

D(cm /sec) 

T(°K) 

934 

1.10*10"q 

1122 

939 

1.41*10 q 

1127 

950 

2.06*10 Q 

1129 

956 

2.20*10“^ 

1133 

966 

2.63*10 Q 

1140 

982 

2.87*10"q 

1145 

990 

5.03*10 q 

1181 

990 

2.63*10 q 

1183 

1002 

6.26*10“^ 

1190 

1008 

4.65*10"q 

1205 

1013 

3.86*10 ^ 

1227 

1022 

6.49*10 q 

1248 

1026 

4,68*10 ^ 

1259 



1272 

1066 

9.99* 10" 


1073 

1.49* IO'q 

1302 

1076 

2.11*10 1 

1323 

1083 

1.71*10"q 

1332 

1087 

1.35*10 : 

1340 

1096 

lAS-lOll 

1344 

1101 

3.56*10 1 

1361 

1111 

4.56*10"^ 

1391 




-9 

-9 

-9 

-9 

-9 

**9 

-8 

-8 

-8 

-8 

-8 

-8 

-8 

-8 


-9 

-9 

-9 

-9 

-9 

-9 

-8 
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TABLE V 

CONCLUDED 

Cs- 

133 


T(;°K) 

2 

D(cm /sec) 
. -10 

T(°K) 

2 

D(cra /sec) 
-9 

T(°K) 

2 

D(cm /sec) 
-9 

T(°K) 

2 

D(cm /sec) 
-8 

917 

8.86*10 r 

1133 

3.50-10 q 

978 

1.83-10 q 

1139 

1.24*10 ° 

941 

1.01*10'“q 

1151 

3.63*10 q 

1010 

2.93*10 q 

1142 

1.46*10"'“ 

954 

l,84*10"q 

1174 

5.99-10“^ 

1022 

3.87*10 q 

1149 

2.08- 10"“ 

973 

4.74*10q 


^ -9 

1028 

8.23*10"^ 

1160 

1.89*10’'“ 

990 

4.05* 10“q 

1220 

2.76*10 q 

1044 

1.08-10 q 

1170 

1.76*10"“ 

1008 

1020 

8.30*10"^ 
7.84*10 o 

1239 

1250 

3.30-lOq 
2.71*10 q 

1050 

9.06*10“^ 

1172 

1186 

1.78*10"“ 

4.39*10"“ 

1037 

1.04*10-« 

1267 

1267 

3.01*10 q 
5.30*10"^ 

1060 

1066 

4.26*10"q 
6.51*10 „ 

1193 

1205 

2.36*10"® 

5.83*10"“ 

1038 

1.52. lO^q 

1284 

7.57*10_^ 

1080 

•mQ 

8.78*10 q 

1222 

6.84*10"® 

1056 

1059 

1.39*10 Q 
2.57*10 q 

1294 

1300 

3.46*10 q 
6.50*10 q 

1083 

1083 

5.20*10"^ 

l.ll*10l« 

1252 

5.22*10"“ 

-9 

1060 

2.42*10*’q 

1326 

5.01*10 q 

1094 

1.16*10 q 

1305 

9.67*10 q 

1087 

2.73*10 q 

1330 

8.22*10 q 

1100 

7.88*10 q 

1311 

8.52*10"^ 

1107 

2.30*10 q 

1332 

9.38*10'^ 

1111 

7.85*10”^ 

1344 

1.77-10:“ 

1107 

2.61*10 q 

1462 

2.20*10"^ 

1116 

1.21-10"“ 

1370 

1.63*10 ° 

1122 

3.60*10 

1511 

6.35*10 

1122 

1136 

1.63*10“^ 

1.86*10 

1395 

2.25*10" 


TABLE VI 

DIFFUSION COEFFICIENTS OF THE ALKALI METALS IN 304 STAINLESS STEEL 




Li- 7 


T(°K) 

2 

D(cm /sec) 
^ «-10 


2 

D(cm /sec) 
-9 

1140 

5.57.10 

1289 

1.30*10 1 

1143 

5.81-10 

1316 

1.55-10 q 

1156 

l.ll-io:^ 

1319 

1.97-10 q 

1175 

1.63-10 1 

1340 

3.78-10 q 

1186 

1.62-10 q 

1344 

4.22-10 q 

1189 

2.06-10 q 

1368 

8.23*10'q 

1220 

3.35-10“^ 

1379 

1427 

6.76-10**^ 

2.71-10'** 


Na-23 


948 

983 

3.85-10'’:J^ 

3.58-io':J: 

1208 

1.69-10 

-9 

1002 

5.29*10:}« 

1227 

1.37*10 q 

1044 

6.11-10 q 

1239 

1.17-10 q 

1106 

1.15-10'^ 

1241 

2.19-lO'q 

1122 

1.09- lO'q 

1272 

4.63-10 q 

1133 

1.19-10 q 

1276 

5.54-10 

1149 

1168 

1.35-10 p 
1.19-10"^ 

1326 

1.07-10"^ 


K-39 


T(°K) 

2 

D(cm /sec) 

T(°K) 

2 

D(cm /sec) 

953 

986 

4.51-lo:“ 
7.47-10 q 

1119 

1134 

1.69-10"q 
1.82-10 q 

1008 

1038 

9.64-10 q 

1172 

1175 

2.26-10": 

2.21-10": 

1040 

1.24*10"q 

1198 

2.56.10": 

1050 

1.34-10 q 

1215 

2.75-10"^ 

1062 

1073 

1.51-10:; 
1.25-10 : 

1305 

5.34-10;“ 
8.73-10 q 
3.06-10": 

1075 

1076 

1.30-10 : 
1.58-10': 

1330 

1389 

1109 

1.36*10': 

1433 

6.73-10"^ 

1116 

2.21-10“^ 




Rb~ 

■31 


1043 

6.28-10”}° 

8.51-10 

5.19-lO'q 

1161 

4.79-10"q 

1060 

1065 

1176 

1182 

4.53-10": 

6.78-10": 

1071 

1.49-10': 

1198 

9.27-10"^ 

1080 

1.27-10 : 


— Q 

1116 

1127 

3.58.10 q 
1.94-10': 

1355 

1404 

1.13-10 : 
3.01-10"^ 
1.11- 10“** 

1133 

2.30.10": 

1445 

1155 

4.46-10" 





TABLE VI CONCLUDED 


Cs-133 


T(°K) 

2 

D(cm /sec) 

Ii°K) 

2 

D(cm /sec) 

960 

1.52*10'^ 

1166 

1.23*10“® 

962 

1.22-10 ^ 

1175 

1.24*10"^ 

975 

2.15-10 q 

1195 

1.39-10"^ 

1016 

2.37*10 1 

1200 

2.10*10"^ 

1038 

1041 

3.76-10 1 
4.03*10 Q 

1205 

2.15-10"^ 

1055 

5.78-10 q 

1287 

3.75-10"^ 

1070 

3.79*10 Q 

1379 

6.43*10 1 

1075 

6.19*10 1 

1381 

7.65*10 1 

1087 

6.79*10 1 

1412 

7.87-10 

1105 

6.09-10"^ 

1427 

1.15*10"“ 

1114 

1.39*10 q 

1431 

9.98*10 1 

1116 

9.42*10”^ 

1439 

8.30*10"^ 

1120 

1.31*10"^ 

1443 

1.00*10"“ 

1136 

1152 

1.16-10“° 
1. 66.10" 

1484 

1.39*10"® 



Table VII 


D and Q for IMPURITY DIFFUSION in 304 STAINLESS STEEL 
o ^ 

928 < T(°K) < 1228 


Impurity D^ Q 


Li-7 

Cl. 86 

± 

2.01) 

X 

10° 

52.96 

± 

0 . 54 

Na-23 

(6.11 

± 

0.56) 

X 

10"^ 

14.22 

± 

0.19 

K-39 

(9.53 

± 

0.67) 

X 

10“^ 

14.09 

± 

0.15 

Rb-8 5 

(3.19 

X 

0.58) 

X 

10"^ 

41.92 

± 

0.38 






-4 




Cs-133 

(6.75 

± 

0.70) 

X 

10 ‘ 

25.04 

± 

0.22 


1228 

< 

T(°K) 

< 

1789 




Li-7 

(1.73 

± 

0.31) 

X 

10^ 

84.42 

± 

0.56 

Na-23 

(1.74 

± 

0.59) 

X 

lo"^ 

74.08 

± 

0.93 

K-39 

(1.72 

± 

0.10) 

X 

10^ 

75.22 

± 

0.18 

Rb-85 

(5.56 

± 

1.75) 

X 

10° 

98.16 

+ 

1.01 

Cs-133 

(5.03 

± 

0.49) 

X 

10“^ 

24.50 

± 

0.26 


2 

D^ is in units of cm /sec 
o 

Q is in units of kcal/mole 
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Table VIII 

D and Q for IMPURITY DIFFUSION in TANTALUM 
o ^ 



990 

< 

T(°K) 

< 

1230 




Impurity 



D 

o 




Q 


Li-7 

(5.15 

± 

0.91) 

X 

10’^ 

20.29 

± 

0,41 

Na-23 

(2.07 

i 

0.27) 

X 

10'^ 

11.90 

± 

0.37 

K-39 

(2.71 

± 

0.50) 

X 

10“^ 

32.19 

± 

0.43 

Rb-85 

(1.22 

± 

0.37) 

X 

10'^ 

36.46 

± 

0.72 

Cs-133 

(8.53 

± 

0.61) 

X 

10‘^ 

15.07 

± 

0.16 


1250 < TC°K) 


Li-7 

(3.38 

± 

0.35) 

Na-23 

(1.13 

± 

0.17) 

K-39 

(1.06 

± 

0.07) 

Rb-85 

(2.33 

± 

0.24) 

Cs-133 

(1.18 

± 

0.12) 


< 1790 


X 

10"^ 

25.49 

± 

0.32 

X 

10-^ 

15.46 

± 

0.47 

X 

10"^ 

22.25 

± 

0,20 

X 

10-' 

14.78 

± 

0.31 

X 

10"^ 

13.48 

± 

0.26 


2 

D^ IS in units of cm /sec 
Q is in units of kcal/mole 
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Table IX 


D and Q for IMPURITY DIFFUSION in ZIRCALOY-2 
o 


900 < TC°K) < 1045 

Impurity Q 


Na-23 

C9.76 

± 

2.99) 

X 

10"^ 

26.44 

± 

0.72 

K-39 

C2.12 

± 

0.31) 

X 

10"^ 

30.92 

± 

0.31 

Rb-85 

(1.11 

± 

0.26) 

X 

10^ 

42.82 

± 

0.55 

Cs-133 

(8.91 

± 

2.96) 

X 

10^ 

52.91 

± 

0.79 


Na-23 

K-39 

Rb-85 

Cs-133 


116.19 ± 3.48 
111.06 ± 0.58 
37.59 ± 0.75 
36.74 ± 0.94 


1345 < TC K) < 1520 


Na-23 

(5.46 

1 

0.75) 

X 

lo^o 

K-39 

(5.35 

± 

0.82) 

X 

10^ 

Rb-85 

(1.11 

± 

0.22) 

X 

10"^ 

Cs-133 

(1.21 

± 

0.30) 

X 

10'^ 


2 

is in units of cm /sec 
Q is in units of kcal/mole 


1045 

< 

T(°K) 

< 

1345 




(1.92 

+ 

0.25) 

X 

10'^ 

40.77 

± 

0.37 

(3.08 

± 

0.41) 

X 

10° 

46.10 

± 

0.38 

(3.53 

± 

0.54) 

X 

10’^ 

25.79 

± 

0.44 

(7.74 

± 

1.01) 

X 

10"^ 

35.03 

± 

0.37 
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both a hexagonal close packed, h.c.p,, and a body centered 
cubic, b.c.c., geometry depending on the temperature, zir- 
caloy-2, which has a composition of “^98. 5? zirconium, will 
probably have a h.c.p. and b.c.c. structure. However, in- 
stead of having a single discontinuity, there is a large 

transition region which possesses its own diffusion char- 

4 3 

acteristics. Pande, Naik, and Agarwala have also observed 
this three region activation energy diagram in their study 
of the diffusion of chromium and iron in zircaloy-2. Their 
transition region held limits of 780 and 980*^C, while this 
work extends the transition region from 772 to 1072°C. This 
large transition region can be explained by a statistical 
average of the diffusivity of both phases present at a par- 
ticular temperature 

When examining the stainless, steels, the results in 
this work did not completely agree with literature values. 

A reason for the discrepancy may involve the fact that the 
composition of stainless steel has a wide range in the per- 
cent abundance of its constituent elements. For example, 

44 

Delong has found that 304 stainless steel changes from a 
face centered cubic, f.c.c., to a body centered cubic, b.c.c., 
structure at 910°C, while the results in this study indicate 
a phase change at 955°C. 

While carbon, deuterium and hydrogen, molybdenum,"^^ 

48 

and sulpher have been studied as they diffused through 
tantalum the diffusion of the alkali metals have not been 
studied in tantalum. Tantalum has a body centered cubic 
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structure, thus there will be no phase transformation. How- 
ever, from the data obtained in this experiment it appears 
that there has been a phenomenological occurence at approxi- 
mately 1000°C. There have been many possible explanations 
suggested, but the one that presents the greatest problem 
in discounting it is that of cold working. For tantalum to 
be processed into the thin ribbons used in this study, tan- 
talum would undergo a great deal of cold working. It is 
possible to remove these defects by annealing, and it so 
happens that the annealing temperature for tantalum is 1000°C. 
The resultant reduction in lattice disorder would tend to 
decrease the diffusion coefficient. However, it would in- 
crease when the temperature was again increased. This is 
exactly what occurs in figures 16 - 20. 

These results provide an insight into the basic dif- 
fusion mechanism in addition to supplying D, D^, and Q. For 
example, the relationship between the diffusion coefficient, 

D, and the ionic radius at a constant temperature can be 
compared, by considering the alkali metal series; Li, Na, K, 

Rb , and Cs, in the same host material. Plots of these five 
impurity ions as observed in this study are shown in Figure 
21 for tantalum, Figure 22 for zircaloy-2, and Figure 23 for 
304 stainless steel. Although there is no phenomenalogical 
relationship, the diffusion coefficient versus mass is some- 
what analoguous to sputtering yield data in which very large 
differences in secondary ion yield vary with atomic mass and 
ionic radius. 
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Each set of figures must be considered separately, as 
they represent completely different crystalline structures 
and alloy compositions. 

Tantalum : Referring to Figure 21 A and B, there is a 

similarity in the pattern of the constant temperature lines. 
Independent of the ambient temperature of the material, the 
diffusion coefficient seems to follow the same pattern of 
maxima and minima with increasing ionic radius. Similar dif- 
fusion patterns can be expected if there were no changed in 
either/or both the diffusion mechanism and the geometric 
structure of the crystal. There has been no reported phase 
change in tantalum, which always possesses body centered 

cubic geometry. Note that as both the ionic radius of the 
o o 

impurity, and T ( K) increase, there appears to be an 

_ 8 2 

upper limit on the order of 10 cm /sec for D. At the same 

time there is a ’’collimating" effect on D as r- goes from 

imp ° 

approximately iX to approximately 1.33A. In an attempt to 
determine the reason for the behavior at these two points, an 
intersititial mechanism for diffusion was assumed and the 
largest sized impurity which would fit through the confining 
lattice ions for a b.c.c. system was calculated. For a b.c.c. 
system the following formula was used 

^impurity ion - ^lattice ion (59) 

where a is the edge size. For tantalum, it was found that 
^imp be less than 1.66 This is extremely close to 
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Figure 21 
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the 1.67 X point on the graph for the Cs^ ion. This implies 

that the lattice structure of the host material plays a role 

in limiting the size of the impurity ion’s diffusion. To 

go beyond this size would require additional energy. The 

o 

region of 1 < < 1.33 shows what might be considered 

inhibited diffusion. This could be caused by the simultane- 
ous diffusion of all the other impurities in this dynamic 
system while being analyzed for just one particular impurity. 

Zircaloy- 2 : A very interesting pattern of diffusion 

can be observed in Figures 22 A, B, and C. Using the 900°K 
line in Fig. 22 A as a reference, in a h.c.p. structure, as 
the ionic radius of the impurity increases, the diffusion co- 
efficient decreases. However, as the temperature is increased 
it is evident that the lattice vibrations cause D to increase 
and the larger impurity ions to diffuse at a greater rate. 

This increase in acceptable ionic radius appears to have an 
upper limit of about 1.47 X which corresponds closely to the 
radius of the Rb^ ion. 

Figure 22 B represents the transition region where 
the crystal structure changes from h.c.p. to b.c.c.. This 
graph is not simply understood as the crystal structure 
changes with temperature. It does appear that the disorder 
induced in the crystal by the phase change increases the 
probability of diffusion. 

Figure 22 C shows the zircaloy-2 system in its b.c.c. 
state. Note that DC1250°K) is less than D(1Q40°K) even though 
the temperature has been increased into the higher temperature 
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phase region. A possible mechanism to qualitatively explain 
this sudden drop in diffusion rate is that as the transition 
region is approached from a lower temperature, lattice defects 

begin to occur before the transition point is actually reached. 
These lattice defects are greatly multiplied in number in 
the phase transition region due to the increased disorder 
in the lattice structure. This increases the impurity dif- 
fusion rate in the phase transition region. However, once 
the phase transition is passed, the crystal lattice became 
more ordered, and many of the lattice defects retreat to 
defect "sinks," such as grain boundaries. This results in 
an initial decrease in the diffusion rate. As the defect 
sinks will release these defects into the crystal lattice 
as the temperature increases, the diffusion rate will again 
begin to increase. 

Figure 22 C shows an interesting occurrence as both 

o o o 

T( K) and r- „(A) increase. In the region r. > 1.47 A, D 

imp ^ ® imp ’ 

is not only relatively invariant as increases, but also 

relatively invariant as the temperature increases. It seems 
that as T increases D is constrained to be of the order of 
magnitude of IQ cm /sec. A partial explanation is proposed: 
if it is assumed that the solute will diffuse by an inter- 
stitial mechanism, a relationship for the largest impurity 
ion accepted before configuration problems are encountered 
is again equation (59) . Using of zirconium, it 

is determined that r^^j^p < 1.47 X. From r^^p = 1.47 X in 
Figure 22C, the diffusion coefficient changes very little 
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with increasing Even increasing the temperature does 

0 

not cause much of an increase in D after r . „ reaches 1.47 A. 

imp 

The changing behavior of D in the range of 1.33 
o 

< < 1.47 could possibly be explained by the fact that 

the zircaloy-2 has a certain percentage of tin, which 

o 

has an ionic radius of 0.93 A. This means that for a given dif- 
fusion jump through four restraining lattice ions, if one is 
a tin ion, impurity ions of larger size will be able to dif- 
fuse through easier than if all four lattice ions were zir- 
conium. This would represent only a minor perturbation as 
the amount of tin is limited. At higher temperatures, this 

o 

represents a negligible effect as gets closer to 1.47 A 

the size plays an increasingly inhibiting effect. 

304 Stainless Steel : Stainless steel 304 presents a 

special challenge for determining diffusion parameters using 
a mass spectrometer. No longer is just a pure metal being 
considered, as in the case with tantalum, or a very dilute 
alloy as with zircaloy-2, but an alloy which has major con- 
stituents of iron, nickel and chromium. The activation energy 
plots have shown the phase change representing a transforma- 
tion from f.c.c. to b.c.c. along with the corresponding 
activation energies. On the other hand. Figures 23 A and B 
represent a more complex picture in regard to D as both, 
and T are varied. 

Figure 23 A represents the f.c.c. structure of the 
stainless steel. As each succeeding temperature line is 
examined, a pattern in the region of approximately 0.70 A 
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and 1.5 A takes place with an apparent change of slope and 
vortex, respectively. This can be explained by the fact that 
as the temperature is increased the lattice vibrations are 
increased so that it is easier for the impurity of a larger 
size to diffuse. There is a limit as to the size that will 
be accommodated. This is primarily controlled by the crystal- 
lographic structure of the host lattice. Stainless steel is 
affected in that there may be unit cells of different host 
ions, representing the major constituents of stainless steel, 
and grain boundaries may differ for specific samples. When 

the effect of Ni in the f.c.c. structure is considered r- 

imp 

should be less than 1.07 X. Configurational limitations do 

not appear to be significant until is approximately 

1.30 X and that the region from l.Q X to 1.3 X is a relatively 

temperature and radius independent region. When r- becomes 

o 

greater than 1.33 A, the larger size ion cannot diffuse as 
readily until more energy is supplied to increase the lattice 
vibrations . 

Figure 23 B represents the b.c.c. structure of 304 

stainless steel where there are again inflection points 

which indicate an area of physical constraints. The point 
o 

'x-l A is caused by the influence of Ni, while the area around 
o 

1.3 A is caused by the influence of Cr (r . should be less 

imp 

, o o 

than 1.23 A) and Fe should be less than 1.34 A). 

After r^^p becomes greater than 1.33 A, D drops unless more 
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energy is given to increase the lattice vibrations. There 
appears to be a significant size effect at approximately 1.67 X. 
The increase in D at lower temperatures is probably due to the 
fact that there are unit cells consisting of different sized 
lattice ions, so that larger ions than normal could diffuse 
through. However, this is a limited effect as T increases, the 
increased lattice vibrations no longer tend to increase D and 
the larger size of the impurity ion becomes a more dominant 
parameter . 

In fact, there are two aspects to note. The first is 

that along with the other samples, as both a-nd T increase 

it appears there may be an upper limit for impurity diffusion 

coefficients in these b.c.c. structures of approximately 
^8 2 

10 cm /sec. The second is that as in the case of zircaloy-2, 
after the phase change the diffusion coefficient drops 
initially and then proceeds to rise with temperature. 

As there are no reported values for the diffusion of the 
alkali metals in 304 stainless steel, it was thought that a com- 
parison between Na-23 and C-14 diffusion in 304 and 316 
stainless steel would be desirable. This would enable one to 
determine whether the results obtained via mass spectrometry 
were reasonable in comparison to the determination of diffusion 
coefficients by other methods for a different element. Table X 
illustrates that the Na-23 diffusion is within an order of 
magnitude of the diffusion of C-14. The Na-23 also seems to 
diffuse at a slower rate than the C-14 as would be expected 
since the Na-23 has a larger ionic radius than C-14. The C-14 
diffusion coefficients were determined by Perkins and Carlson, 
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ainless Steel 


Table X 




^-Lxrusion in 304 st 


C *14 Diffusion 


873 

923 

943 


4. 30' 


£Cref , 

^•27*10“^° 
1 .66*10”^® 
3 . 02 


963 

^73 6,54.j^q-10 


^•57'1o“^® 


T 


^Na-23 Diffusion 

^CThis Stilly) 


^^8 3 . 85 . 10 "^° 


993 


8.71-10"^^ 

983 

3.58*10"^° 

1023 


1.91.10"^ 

1002 

3.29-10"^^ 

1049 

1063 

1093 

3.36*10'^ 

3.98.10'^ 

7.25-10"^ 

1044 

3.11*10"^^ 

1123 

9.54*10'^ 

1 . 02 - 10 "® 

1106 

1122 

i.15'10"^ 

l-09-lo"^ 

1153 


2 . 19 . 10 "® 

1133 

1149 

i-l9-l0“? 

1-35-10"^ 

1198 


^.36*10"® 

1168 

1.19-10"^ 

1209 2 

1223 

.83*10"® 

6 . 61 . 10 "® 

1208 

1.69-10"^ 

1269 9 

1273 

74'10"® 

i.45-lo"7 

1227 

1239 

1241 

i-37-10"^ 

l-17-10"f 

2.19-10"^ 

1333 

1373 

1473 


2. 75 . 10" 7 
5 . 01 - 10 "' 
1 . 54 - 10 "® 

1276 

1326 

5.54-l0"f 

1.07-10"® 
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and by Agarwala, Naik, Anand, and Paul.^^ 

In the further analysis of the data, two additional 
relationships merit examination. Based on equation (48) 

D = exp (-Q/RT) 

it would seem that there should be some type of direct re- 
lationship between and Q. Using the PDF 8/1 computer a 
best fit least squares program was run utilizing data on 
and Q for h.c.p., f.c.c., and b.c.c. structures and at all 
temperatures. It was found that the correlation coefficient 
for the best fit line was 0.98 where 1.0 is a perfect fit. 
Subsequently, it was determined that the relation between 
and Q was 

= 3.3 X 10'^ exp {(0.36)(Q)} . (60) 

Whether or not this equation applies only to a dilute solute 
diffusion of the alkali metals in these particular matrices, 
or for any other diffusing system requires further investiga- 
tion . 

A second interesting consequence of this work involves 
the possibility of using this method to further develop 
concepts on ultra-fast diffusion. Ultra-fast diffusion 
has been considered when the solute diffusion is much greater 
than self -diffusion in the host material. An example is the 
diffusion analysis of zircaloy-2 in this work. The self-dif- 
fusion rates of pure a-zirconium have been determined by 
Dyment and Libanati.^^ For a temperature range of 1036°K to 

1130°K, they found D to be on the order of 10~^^ to 10~^^ 

2 

cm /sec, while in this study, the diffusion rates of the alkali 
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metals for the same temperature range were on the order of 
10'^ to lO”^ cm^/sec. This is approximately four to five 
orders of magnitude faster. Ultra-fast diffusion has been 
suggested by Hood and Schultz to take place by an inter- 
stitial mechanism. They studied both Co-60 and Fe-59 dif- 
fusion in pure a- zirconium and found that D was of the order 
of 10' to 10 cm /sec which corresponds well with this 
work. That ultra-fast diffusion would take place by an 
interstitial mechanism coincides with the mechanism for de- 
termining the largest diffusing ion size; a concept used to 
explain figures 22 and 25 . 

It should be further noted that in Part II, con- 
cepts of the jump frequency w , the mean vibrational frequency, 
V , and the change in free energy, aG, were all discussed. 

With the results obtained in this study, it now becomes im- 
portant to determine some of these quantities. 

59 

Zener has stated that while it would not be 
impossible for aG to increase with increasing temperature, 
it would be highly improbable. The reason for this is that 
a negative temperature coefficient of the elastic moduli 
exists. Zener also found that the change in entropy, AS 
can be related to the activation energy, Q, by 

AS - ABQ/T^ (61) 

where x is a aiumerical coefficient essentially equal to 
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unity for interstitial, 6 is a dimensionless quantity tabu- 
lated by Zener for some pure metals, and is the melting 
temperature for the host material. The value of AS can be 
used to calculate v using 

Dq = (1/6) V exp (4S/R) (62) 

provided that is known. 

Of the three sample materials used in this study, 

6 has been determined for tantalum only. It has a value of 
0.40. Using values for obtained experimentally, v can be 
determined. Once the values of v have been determined, co 
can be calculated using equation (23) . 

0 ) = vN^ = V exp (-Q/RT) exp (aS/R) 

Table XI gives the calculated values for u for the five 
diffusing impurity isotopes in both temperature regimes. 
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Table XI 


Calculated Values o£ v and to For Tantalum 


990 < T(°K) < 1230 


Impurity 

V 

(sec’^) 

exp C-Q/RT) 

exp CaS/R) 

0 ) (sec 


Li-7 

2.83 

X 


.997 

1.001 

2.82 

X 

10^^ 

Na-23 

1.14 

X 


.998 

1.0007 

1.13 

X 

10^® 

K-39 

1.49 

X 


.995 

1.002 

1.48 

X 

10^^ 

Rb-85 

6.71 

X 


.994 

1 . 002 

6.68 

X 

10^^ 

Cs-133 

4.69 

X 


.998 

1.001 

4.68 

X 

10^® 




1230 

< T(°K) < 

1790 




Li-7 

1.86 

X 

10^^ 

.996 

1.002 

1.85 

X 

10^^ 

Na-23 

6.22 

X 

10® 

.998 

1.001 

6.21 

X 

10® 

K-39 

5.83 

X 

10^® 

. 997 

1.001 

5.82 

X 

10^® 

Rb-85 

1.28 

X 

10® 

.998 

1.001 

1.27 

X 

10® 

Cs-133 

6.50 

X 

10® 

.998 

1.001 

6.49 

X 

10® 
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Part VI 


CONCLUSIONS 

It was the purpose of this work to illustrate a new 
technique in the general field of materials characterization. 
The transport of impurities, in the p.p.m. to p.p.b. range 
and less, is becoming increasingly important. Since the 
mass spectrometer is an instrument which can detect impurity 
concentrations on this scale, it offers another view of the 
high temperature diffusion properties in materials which are 
of engineering usefulness. This is instead of being confined 
to the world of the single crystal. 

As a consequence of this work, not only were the dif- 
fusion coefficients of the alkali metals in the sample 
materials calculated, but known phase transformations were 
determined in the samples of 304 stainless steel and zircaloy- 
2. This, of course, lends a measure of credibility to the 
results obtained. 

While the primary intent of this work has been accom- 
plished, it has been shown that what used to be a chemists 
tool can be applied in metallurgical research. It is not 
to be suggested that the mass spectrometer is the ultimate 
instrument, but rather that it allows one more unique way of 
looking at material's properties on the atomistic level. It 
is conceivable that by determining the relationship between 
the diffusion coefficient and the atomic radius of all the 
impurities in a material that a better understanding of 
diffusion processes could be obtained. It is hoped that 
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further funding will allow this laboratory to pursue this 
and to see how ion implantation of selected species affect 
the diffusion process in certain materials - i.e., semi- 
conductors . 
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I 


APPENDIX 1 


Program to Determine Least Squares Fit o£ 

Data to Exponential Curve, Correlation Coefficient, 
and 95% Confidence Limit 


C-8K FOCAL ^11969 

01.10 T PROGRAfv) FITS PAIRS OF POINTS TO AN EXPONEiNT I AL 

01.11 T !"LEAST SQUARES CURVE FIT - TO THE EQUATION: Y=A*£XP( ^ ) 

01.12 T X IS The rectangular absyssa 

01.13 T !" Y IS THE RECTANGULAR ORDINATE (INTERNAL CONVERSION TO LOCI' 

01.14 T !** A IS THE Y-INTERCEPT 

01.15 T !'* B IS THE SLOPE ON SEMI-LOG PAPER 

01.16 T !•’ R IS THE CORRELATION COEFFICIENT 

01.17 T !'• CL IS THE 9 5% CONFIDENCE LIMIT +/- AT X=(TEST POINT) 

01.18 T ! !''PR0GRAM operation: 

01.19 T ! *■ AFTER FIRST PROMPT ENTER THE NUMBER OF ORDERED PAIRS 

01.20 I !■• then -enter X S y AS REQUESTED BY THE PROGRAM 

01.21 T !” -enter rC0.vl25) FOR THE SPECIFIED D.F. vjHEN REQUESTED 

01.22 T !•’ -ENTER THE NUMBER OF TEST POINTS 

01.23 T !" -ENTER THE TEST POINTS 

01.30 A !!!''EMTER THE NUMBER OF VALUES: ”N 

01.35 S SX = 0;S SY = 0J$ XT = 0;S XY=:0;S YT=0 
01 .40 I (N-1 ) 9.9> 1 .6 

01.50 F 1= I jN; D 2 

01.55 G 3.1 

01.60 T !"TRV again -- 1 IS NO GOOD 

01.65 G 1.3 

02.10 A !"ENTER X: "X, "ENTER Y: ",Y;S Y=FLOG(Y) 

02.15 S SX = Sx + xjS SY = SY + YjS XT = XT + Xt2;S XY=XY+X4«Y;S YT=YT+Yt2 

03.10 S XS=SX/N;S YS=SY/N 

03.15 T !!"THE AVERAGE VALUE OF X IS: ",%,X$ 

03.20 S YY=YT- (Syt2/N ) 

03.30 S XX = XT- ( SX t2/N ) 

G3.40 S TT=X Y- ( SY*SX/N> 

03.50 S SS=(TTt2)/XX 

03.55 S MS=(YY-SS)/(N-2) 

03.60 C CALCULATION OF THE SLOPE -- □ 

03.65 S B =T T / X X 

03.70 C CALCULATION OF THE Y-INTERCEPT -- A 
03.75 S A=FExP(VS-B*X$) 

03.80 C CALCULATION OF THE CORRELATION COEFFICIENT -- R 
03.85 S RsTT/FSOT C XX+YY ) 

03.90 G 5.10 
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0-6.10 C PRINT OUT THE RESULTS 

04.20 T %,!!! "RESULTS : 

04.30 I !" THE Slope is: ">b 

04.40 r !" THE r-INTERCEPT IS: ",A 

04.50 T !," THE CORRELATION COEFFICIENT IS: ",R 

04.60 R 

04.70 T !!!!"IF YOU WISH TO STOP -- AT THE NEXT PROMPT TYPE <-l) 

04.80 T !" OTHERWISE CONTINUE 

G4 • 9 0 G 1.3 

05.10 C INPUT FOR THE DE T ERM I N AT I ON OF THE 95% C.L.’S FOR EACH TEST PT 

05.15 T T.2.0> ' .'"ENTER THE VALUE OF 7(0.025) FOR ",N-2, " DEGREES 

05.16 A " OF FREEDOM: ",DF 

05.20 A !"ENTER THE NUMBER OF TEST POINTS: "jM 
05.25 F 1=1, m;D 6 

05.30 DO A 

05.35 F 1=1, M;S CL ( I ) =DF *FSG T C MS* ( 1 /N+ ( x ( I ) - 0 ) t 2 /XX ) ) 

05.36 F 1=1, M;S Y ( I )=FEXP(B*X( I > )*A 

05.37 F 1=1 ,M;D 7 
05.35 0 9.04 

05.39 D 9.05 

05.40 F 1=1 ,M?0 9.1 
05.45 G 4.70 

06.05 C INPUT OF TEST POINTS -- X(/vi) 

06.10 T %2.0,!"ENTER TEST POINT X(",I,"): " 

06.20 A X ( I ) 

07.10 S ZA = FEXP(FLOGC Y( 1 ) ) +CL( I ) ) 

07.20 S ZB=FEXP(FL0G( Y( I ) )-CL( I ) ) 

07.30 S ZC( I ) =FAeS ( ZA-ZB )/2 

09.04 T .'.'"TEST POINTS AND THEIR RESPECTIVE 95% C-L-'S 

09.05 T ! !" M ","X(M> Y(M) CLtI 

09.10 T .' ,%2 .0> I ",X(I)>" ",Y(I)," ",ZCCI) 

09.90 T !"TH£ END"jQ 
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APPENDIX 2 


Program to Determine D for a Given Temperature 
Provided D^ and Q are Supplied 


C-8K FOCAL 9 1969 

01.01 C PROGRAM: DIFF 

01.02 C 

01.10 A !IJ"WHAT IS THE SAMPLE MATERIAL? *'#A 

01.20 A !"WHAT is the diffusing ISOTOPE? ’*,A^B 

01.30 A !*'WHAT IS THE TEMPERATURE RANGE? '*,A,'*TO B > ’‘DEGREES KELVIN." 

01.35 A '"WHAT IS THE VALUE OF D(0)? "DC 

01. A0 A !"WHAT IS THE VALUE OF THE ACTIVATION ENERGY? "Q 

01.60 T 

02.10 A !"T= ">T 
02.15 I (T> 2. 6,2. 5 

02.20 S D=D0^FEXPC-Q/(2E-3)*T) 

02.25 G 3.1 

02.30 T " ","D= "D 

02.31 G 2.1 

02.35 T " ","D= "D 

02.40 G 2.1 

02.50 G 1.1 

02.60 T ! J ! ! ; 0 

03.10 I C 1000-T)2.3>2.3,2.35 
* 
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